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ABSTRACT 

In  this  work,  we  look  at  single  user  and  multiuser  Multiple-Input 
Multiple-Output  (MIMO)  beamforming  networks  with  Channel  Dis¬ 
tribution  Information  (CDI).  CDI  does  not  need  to  be  updated  every 
time  the  channel  changes,  but  only  when  the  statistics  of  the  channel 
change.  Thus,  feedback  in  the  network  is  significantly  reduced  when 
compared  to  CSI  schemes.  With  statistical  information,  we  can  only 
guarantee  quality  of  service  for  a  certain  probability  of  outage  in  the 
network.  The  optimal  beamformers  and  an  optimal  power  control 
algorithm  to  minimize  the  power  cost  function  in  the  network  are 
presented  for  the  well-known  Kronecker  model. 

Index  Terms —  covariance  feedback,  beamformers,  multiuser 
MIMO,  outage  probability,  Channel  Distribution  Information 

1.  INTRODUCTION 

Multiuser  multiple-input  multiple-output  (MU-MIMO)  networks 
have  generated  much  research  interest  in  recent  years.  Ideally,  to 
achieve  maximum  throughput  and  reliability,  all  the  nodes  need  to 
have  perfect  Channel  State  Information  (CSI)  of  all  the  links  in 
the  MU-MIMO  network.  With  this  knowledge,  utilizing  intelligent 
beamforming,  interference  in  the  system  can  be  minimized  and 
optimal  rates  can  be  achieved. 

This  assumption  may  be  feasible  in  small  networks  where  the 
channels  change  slowly.  However,  channels  change  constantly  in 
mobile  networks.  Thus  much  work  has  focused  on  how  to  achieve 
good  system  performance  with  reduced  feedback.  One  transmission 
scheme  uses  limited  CSI,  where  nodes  feed  back  their  CSI  only  to 
nearby  nodes  [1].  This  approach  assumes  that  the  nearby  nodes  dom¬ 
inate  the  interference,  and  should  therefore  yield  near-optimal  per¬ 
formance.  However,  even  in  point-to-point  systems,  low  channel 
coherence  times  make  full  CSI  feedback  infeasible.  Many  schemes 
thus  utilize  quantized  forms  of  CSI  to  reduce  feedback  [2].  However, 
these  approaches  still  suffer  from  having  to  feed  back  information 
every  time  the  channel  changes. 

An  alternate  approach  is  to  utilize  the  channel  statistics,  or 
Channel  Distribution  Information  (CDI),  to  enhance  communica¬ 
tion.  Because  it  takes  into  account  the  randomness  in  the  channel, 
CDI  is  more  robust  to  small  channel  coherence  times  and  is  thus 
valid  for  much  longer  than  CSI.  In  addition,  statistical  data  based  on 
node  location  can  also  be  collected  and  stored  a  priori,  eliminating 
the  need  for  real-time  channel  feedback.  The  ergodic  capacity- 
optimal  input  covariance  for  the  SU-MIMO  channel  with  mean  and 
covariance  channel  information  is  presented  in  [3].  This  result  is 
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extended  to  the  MIMO  multiple  access  channel  in  [4].  Suboptimal 
solutions  for  the  MIMO  broadcast  channel  are  given  in  [5].  These 
works  focus  on  average  capacity,  but  in  doing  so,  do  not  account 
for  a  fundamental  issue  with  using  CDI:  much  of  the  time  the 
throughput  will  be  significantly  less  than  the  average  due  to  varying 
channel  conditions.  Thus,  previous  schemes  do  not  provide  robust 
transmission  over  time. 

This  work  seeks  to  address  this  problem  by  looking  at  the  outage 
of  the  links  in  the  network.  Due  to  the  randomness  in  the  channel, 
reliable  transmission  cannot  be  achieved  all  the  time  using  just  CDI. 
However,  in  this  work  we  guarantee  a  certain  signal-to-interference- 
plus-noise  ratio  (SINK),  and  thus  reliable  transmission,  for  a  speci¬ 
fied  probability  on  all  the  links  using  power  control  and  linear  beam¬ 
forming.  While  many  works  in  the  literature  assume  the  channels  ex¬ 
perience  independent  identically-distributed  (i.i.d.)  Rayleigh  fading 
for  analytic  convenience,  the  Kronecker  model  captures  correlation 
between  the  antennas  and  is  a  more  accurate  representation  of  what 
happens  in  real  systems  [6].  With  the  assumption  of  the  Kronecker 
model  on  all  the  links,  the  expression  for  the  outage  probability  on 
each  link  is  derived.  Then,  we  derive  the  optimal  transmit  and  re¬ 
ceive  beamformers  given  this  model  and  give  an  optimal  power  con¬ 
trol  algorithm.  In  Section  2,  the  problem  formulation  is  given,  and 
the  expression  for  outage  probability  is  derived.  Section  3  derives 
the  optimal  beamformers  and  gives  the  power  control  algorithm.  We 
plot  our  derived  CDI  solution  for  different  outage  thresholds  in  Sec¬ 
tion  4,  and  we  summarize  our  results  in  Section  5. 

In  this  work,  the  following  notation  is  used:  italicized  letters  in¬ 
dicate  scalars  (e.g.  pi),  lower-case  bold  letters  indicate  vectors  (e.g. 
v),  and  upper-case  bold  letters  indicate  matrices  (e.g.  A).  Further¬ 
more,  indicates  the  Hermitian  operator,  and  (•)*  indicates  the 
conjugate  transpose  operator. 

2.  PROBLEM  FORMULATION 
2.1.  System  Model 

This  work  considers  time-varying  MIMO  channels  for  many  users 
in  a  network.  Consider  a  MIMO  network  with  L  transmit-receive 
pairs.  At  link  I,  the  transmitter  sends  the  symbol  xi  (t)  to  the  receiver. 
The  transmitter  uses  unit-norm  beamforming  vector  v;  (t)  to  precode 
the  signal,  and  transmits  with  power  pi  {t) .  The  receiver  employs 
the  linear  unit-norm  beamformer  u;  (t)  to  combine  the  signal.  The 
channel  from  transmitter  i  to  receiver  I  is  given  by  Hud).  The  noise 
Ni  (t)  is  distributed  as  a  complex  circular  Gaussian,  and  represents 
the  combined  noise  after  applying  the  receive  beamforming  vector 
to  the  incoming  signal. 

In  schemes  that  use  perfect  CSI,  for  every  change  in 
all  the  transmit  and  receive  beamformers  must  be  updated,  and  the 


power  allocation  scheme  changes.  This  work  will  consider  the  case 
where  the  channel  varies,  but  the  statistics  stay  constant.  There¬ 
fore  the  power  allocation  and  transmit  and  receive  beamformers  stay 
fixed,  and  the  time  variable  will  be  dropped  for  notational  conve¬ 
nience.  To  further  simplify  notation,  define  Gh  =  luf^H^Vip  as 
the  beamforming  channel  gain  from  the  transmitter  on  link  i  to  the 
receiver  at  link  I  and  as  the  noise  power  for  the  Z’th  link.  Then, 
under  this  model,  the  SINR  T i  on  each  link  can  be  shown  to  be 


r. 


PiGii _ 


(1) 


If  perfect  CSI  is  available,  to  ensure  a  reliable  link  is  available 
to  all  nodes  in  the  network,  each  link  has  an  SINR  constraint:  T; 
must  be  greater  than  a  threshold  7; .  The  goal  is  then  to  minimize  the 
power  consumed  by  the  network  while  meeting  all  the  SINR  con¬ 
straints  [I].  When  only  CDI  is  available,  however,  the  Hk’s  are 
unknown.  Instead,  they  are  assumed  to  be  a  random  variable  drawn 
from  a  complex-normal  distribution,  with  a  Kronecker  structure  on 
their  correlation.  The  Kronecker  model  assumes  the  transmit  and  re¬ 
ceive  correlations  are  independent  at  the  link  ends.  This  implies  that 
each  transmitter  has  its  own  transmit  correlation  matrix,  regardless 
of  the  receiver,  and  each  receiver  has  its  own  receive  correlation  ma¬ 
trix,  regardless  of  the  transmitter.  Thus,  the  channel  from  transmitter 
i  to  receiver  I  is  given  as: 


Proof:  First,  consider  Z  =  u^Hv.  Z  is  a  linear  combination  of 
0-mean  complex  Gaussians,  so  Z  is  0-mean  complex  circular  Gaus¬ 
sian.  The  next  step  is  to  find  the  variance.  To  do  this,  perform  the 
following  calculations: 


var{Z)  =  E[Z  ■  Z*]  =  £;[(u"Hv)(u^Hv)*] 

=  i!l[vec(u^Hv)vec(u^Hv)^] 

=  E[{v'^  0  u‘^)vec(H)vec(H)‘^(v*  0  u)] 

=  (v^STv)(u^SfiU) 

This  gives  an  expression  for  var{Z).  The  random  variable  of 
interest,  however,  is  \Z\‘^.  The  norm  squared  of  a  unit  variance  com¬ 
plex  circular  Gaussian  has  a  7I,  or  exponential,  distribution.  Then 
jZj^  needs  to  be  normalized  appropriately,  which  can  be  done  by 
dividing  by  the  variance,  giving  the  result. 


□ 


To  simplify  notation,  define  f  —  (v[^StiVi)  and  r;  = 
(uPSij, u;),  I  —  I, . . . ,  L.  Without  loss  of  generality,  the  ex¬ 
pression  for  outage  probability  of  the  first  user  will  be  shown.  The 
theorem  can  then  be  stated  as  follows: 


vec(H,,)  ~C'Ar(0,S?. 

The  knowledge  of  the  transmit  and  receive  correlation  matrices  com¬ 
prise  the  CDI  of  the  network.  Under  this  model,  the  SINR  becomes 
a  random  variable  since  it  depends  on  the  channel.  Therefore,  the 
constraints  become  outage  constraints,  so  links  are  allowed  to  have 
an  SINR  below  their  thresholds  7;  for  some  probability  a; .  The  cost 
function  considered  in  this  work  is  the  weighted  sum  power.  In  this 
setup,  each  link  I  in  the  network  incurs  some  cost  wi  >  0  to  transmit 
across  its  link. 

To  simplify  notation,  define  the  weighting  and  power  vectors 
as  w  =  {wi, . . .  ,wl}  and  p  =  {pi, . . .  ,Pl},  respectively.  The 
beamforming  matrices  are  defined  as  U  =  {ui, . . . ,  u_t}  and  V  = 
{vi, . . . ,  Vi}.  The  main  optimization  problem  for  using  CDI  can 
then  be  formulated: 

min  w^p  (2) 

p>0,U.V 

s.t.  Pr(r;  <  7;)  <  a;,/ =  1, . . . ,  T 

To  solve  this  problem,  a  closed-form  expression  for  the  outage  prob¬ 
ability,  Pr(r;  <  7;),  must  first  be  derived. 

2.2.  Derivation  of  Outage  Probability 

In  order  to  get  a  closed-form  expression  for  the  outage  probability, 
r ;  is  first  expressed  as  a  ratio  of  an  exponential  random  variable  to  a 
weighted  sum  of  exponential  random  variables  plus  some  constant. 
The  CDF  for  F ;  is  then  given  in  [7],  which  gives  the  final  expression 
for  the  outage  probability. 

First,  the  following  Lemma  is  needed: 

Lemma  1  Ifu  and  v  are  unit-norm  vectors  and  H  is  Gaussian  ma¬ 
trix  distributed  as  vec(H)  ~  CN[Q,  Sy  (g)  Sh),  then 

_ _ (3) 

(v"Stv)(u^^Sku)  ^ 


Theorem  1  In  a  MU-MIMO  network  where  all  links  experience  cor¬ 
related  Rayleigh  fading  under  the  Kronecker  model  and  the  transmit¬ 
ter  and  receiver  both  employ  linear  beamforming,  the  expression  for 
outage  probability  for  the  SINR  is  given  by 

2 

pout  =  Pr(Fi  <  71)  =  1  - 

Proof:  To  apply  the  results  from  [7],  pout  needs  to  be  in  the  follow¬ 
ing  form: 

Pout  =  Pr(ri  <  71)  =  Pr(^^^^  <  71),  (5) 

where  X  is  an  exponential  random  variable  and  U  is  a  weighted  sum 
of  independent  exponential  random  variables.  Applying  Lemma  I  to 
the  Gu  terms  in  the  expression  for  Fi  in  Eqn.  (1)  gives 

_  riGpiZi 

Ef=2  ’ 

Zi  ~  xi^'si.i.d. 

Now,  divide  the  top  and  bottom  of  the  right  hand  side  of  the  equation 

by  rifipi  and  define  ki  =  Iff}’  ^  ~  Z\,  Y  =  ^1^2  ^iZi,  and 
^2 

(7^  =  .  After  plugging  in  these  variables  into  the  expression 

for  Fi,  pout  has  exactly  the  desired  form  of  Eqn.  (5).  The  closed- 
form  expression  given  in  [7]  can  then  be  applied  to  (5),  which  gives 
the  final  result  for  pout  above. 

□ 

To  apply  this  result  to  (2),  first  define  Pi  =  (1  —  0;;)“^.  Then, 
manipulating  the  constraints  gives  the  following  optimization  prob¬ 
lem: 


s.t. 


T 

mm  w  p 

p>0,U,V 


21^ 


<  /3i,  Z  =  1, . . . ,  L 


(6) 


This  expression  can  also  be  generalized  to  channels  with  general 
correlation  structures  (without  the  Kronecker  model  restriction),  and 
is  derived  and  discussed  in  [8]. 


the  optimal  solution  [9].  Then,  to  find  the  optimal  p,  first  take  the 
logarithm  of  the  constraint  equations  in  (6)  and  multiply  both  sides 
by  ■  This  results  in 


2riti  log  pi 


+ 


Pi 

log  Pi 


^log 


<  Pi 


Now  define 


3.  DERIVATION  OF  THE  OPTIMAL  SOLUTION 

3.1.  Optimal  Transmit  and  Receive  Beamforming 

First,  focus  on  the  optimal  transmit  beamformers.  Define  qi  — 
tiPiyl  =  By  noting  that  p;  =  the  following  equiv- 

alent  optimization  problem  to  (6)  is  presented: 


L 

Ewiqi 

— 

—  ■  '  1=1  ‘ 

s.t.  e^’-iii  <  Pi,l=l,...,L 

V  9!/ 


(7) 


In  the  above  optimization  problem,  the  ti ’s  only  appear  in  the  opti¬ 
mization  function,  and  not  in  the  constraints.  To  minimize  the  objec¬ 
tive  function  in  (7),  then,  the  ti  terms  should  be  maximized.  Since 
ti  is  a  function  of  v; ,  maximizing  ti  is  equivalent  to  solving  the  fol¬ 
lowing  optimization  problem: 


maxvf^STjV;  (8) 

s.t.  ||v;||2  =  l 


7i(p) 

I(P) 


2riti  log  Pi  log  Pi  ^ 
[7i(p),...,/l(p)] 


Upi 

tiPi 


The  function  I(p)  is  identical  to  the  function  presented  in  [9],  with 
different  constants,  and  in  their  work  is  shown  to  be  a  standard  inter¬ 
ference  function.  A  key  property  of  standard  interference  functions 
is  that  they  satisfy  p  >  I(p).  Thus,  using  this  function  and  a  starting 
power  vector  p,  the  update  equation  for  the  algorithm  is  given  as: 


p(-+i)  =  (9) 

The  complete  algorithm  can  then  summarized  as  follows: 

Algorithm  1:  Joint  Power  Control  and  Beamforming 

1 .  Solve  for  the  optimal  u; ’s,  and  vi’s,l  —  1 , . . . ,  L,  as  the  nor¬ 
malized  dominant  eigenvectors  of  Sij, ’s  and  respec¬ 

tively,  as  shown  in  Section  3.1 

2.  Initialize  p  >  0  and  update  p  using  (9)  until  convergence 

In  the  special  case  of  SU-MIMO  (where  L  =  1),  the  optimal 
beamformers  have  the  same  solution  as  in  the  MU-MIMO  case,  and 
the  optimal  power  p  can  be  given  in  closed  form: 


Popt  — 


_ _ 

2(v"tSTVopt)(u"tSflUopt)  log(^) 


(10) 


This  optimization  problem  has  a  well-known  solution:  v;  is  the  nor¬ 
malized  eigenvector  corresponding  to  the  maximum  eigenvalue  of 
St; ,l  =  1,  ■  ■  ■  ,L.  Thus,  this  gives  the  optimal  V  to  solve  (7). 
Since  the  problem  in  (7)  is  equivalent  to  the  problem  in  (6),  the  op¬ 
timal  transmit  beamformers  V  will  be  the  same  as  well. 

Now  focus  on  the  optimal  receive  beamformers.  In  (6),  note  that 

ri  only  appears  in  the  exponential  term,  e  7  p; ,  of  the  constraint 
equation.  The  objective  function  will  be  minimized  when  the  pi ’s  are 
as  small  as  possible  while  still  meeting  the  constraints.  In  order  for 
Pi  to  be  minimized  while  meeting  the  constraint,  observe  that  the 
exponential  term  with  respect  to  r;  should  be  minimized.  Since  r; 
appears  in  the  denominator  of  the  exponential  term,  this  term  will  be 
minimized  when  r;  is  maximized.  Thus,  the  optimal  receive  beam- 
former  for  the  P^  constraint  will  be  the  u;  that  maximizes  rj.  This 
is  analogous  to  the  optimal  transmit  beamformer  case,  so  ui  is  the 
normalized  eigenvector  corresponding  to  the  maximum  eigenvalue 
of  ,l  —  1, . . .  ,L.  This  gives  the  optimal  receive  beamformers 
U. 

3.2.  Power  Control 

In  the  following  analysis,  the  transmit  and  receive  beamformers  are 
assumed  to  be  fixed.  It  can  be  shown  that  if  all  the  beamformers 
are  fixed,  the  constraints  in  (6)  are  satisfied  with  equality  and  thus 
uniquely  determine  the  optimal  power  vector  for  all  choices  of  w  in 


4.  RESULTS 

We  conduct  some  numerical  experiments  to  understand  the  efficacy 
of  the  joint  power  control  and  beamforming  algorithm  developed 
above.  Results  are  given  for  a  single-user  MIMO  setup  and  a  MU- 
MIMO  setup.  All  transmitters  and  receivers  have  4  antennas,  and 
the  SNR  is  held  constant  at  10  dB.  Equal  SINR  thresholds  and  equal 
outage  constraints  are  considered  for  all  users  in  the  system. 

For  single-user  MIMO,  an  algorithm  using  perfect  channel 
knowledge  is  compared  to  the  algorithm  presented  here,  which  uses 
the  covariance  information  only,  in  Experiment  1.  The  covariance 
matrix  here  uses  the  Kronecker  model  parameters  extracted  from 
a  covariance  matrix  generated  from  an  angular  spread  model  using 
100  scatterers  with  a  transmit  and  receive  angular  spread  of  ^  com¬ 
ing  in  at  broadside.  The  algorithm  using  perfect  channel  knowledge 
uses  the  principal  left  and  right  singular  vectors  for  receiver  and 
transmitter  beamforming,  respectively,  and  calculates  the  minimum 
transmit  power  required  to  achieve  the  threshold  (this  is  an  optimal 
scheme).  See  Fig.  1. 

The  plot  yields  some  interesting  results:  the  average  transmit 
powers  for  the  covariance  schemes  are  actually  lower  than  for  using 
the  true  channel  at  5%,  10%,  and  20%  outage  values.  To  calculate 
this,  only  the  covariance  information  is  required,  so  it  needs  to  be 
updated  only  when  the  channel  statistics  change.  Note  also  that  the 
plots  are  in  log  scale,  so  the  savings  over  true  channel  knowledge 
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Fig.  1.  Single  User  Case 
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Fig.  2.  Two-User  Case 


are  significant;  at  5%  outage,  there  is  roughly  a  3dB  gap,  and  at 
10%  outage  there  is  roughly  a  6dB  gap.  The  tradeoff  is  that  the 
covariance  algorithm  will  allow  the  link  to  fail  with  the  specified 
outage  probability.  However,  this  may  be  a  desirable  tradeoff  for 
many  applications  since  it  prevents  the  transmitter  from  attempting 
to  increase  the  power  to  compensate  for  poor  channel  conditions, 
thereby  prolonging  battery  life. 

For  the  MU-MIMO  setup,  a  network  with  2  transmit/receive 
links  and  a  weighting  vector  w  =  [10, 1]^  is  considered.  The  covari¬ 
ance  matrices  use  the  Kronecker  model  parameters  extracted  from  a 
covariance  matrix  generated  from  an  angular  spread  model  using  100 
scatterers  with  a  transmit  angular  spread  of  ^  and  a  receive  angular 
spread  of  For  link  1,  the  beams  are  centered  at  broadside,  and 
for  link  2,  the  beams  are  centered  at  the  null.  In  Experiment  2,  dif¬ 
ferent  outage  thresholds  obtained  from  Algorithm  1  are  compared. 
See  Fig.  2. 

As  expected,  this  experiment  shows  that  as  the  outage  probabil¬ 
ity  threshold  is  increased,  the  power  required  by  the  network  to  meet 
the  threshold  also  increases.  Also,  for  the  multiuser  case,  the  solu¬ 
tions  generated  by  Algorithm  1  become  infeasible  (when  the  curves 
go  to  oo)  as  SINK  requirement  increases.  The  lower  the  outage,  the 
lower  the  SINR  threshold  has  to  be  for  the  outage  to  become  infea¬ 
sible.  The  SINR  thresholds  here  are  low  due  to  the  limitations  of 
the  Kronecker  model  and  the  inability  of  the  optimal  beamformers 
to  suppress  interference,  as  discussed  in  Section  3. 

5.  CONCLUSION 

In  this  work,  a  framework  has  been  presented  for  analyzing  MU- 
MIMO  networks  when  only  the  Kronecker  model  parameters  are 
known.  The  expression  for  outage  probability  was  derived  in  closed- 
form,  and  the  optimal  solution  for  the  beamformers  and  power  con¬ 
trol  was  given.  When  compared  to  CSI,  the  computational  complex¬ 
ity  for  the  algorithm  and  feedback  information  required  for  covari¬ 
ance  information  is  drastically  reduced  since  both  must  be  done  only 
once  for  valid  covariance  information.  Also,  in  the  single-user  case. 


the  amount  of  power  used  at  the  receiver  is  reduced  when  compared 
to  the  optimal  algorithm  using  CSI  on  average. 

The  assumption  of  the  Kronecker  model  in  the  multiuser  case  is 
not  realistic  in  many  scenarios.  The  scatterers  between  one  receiver 
and  transmitter  is  likely  to  differ  from  another  receiver  and  the  same 
transmitter,  and  so  the  correlation  structures  can  differ  at  both  the  re¬ 
ceive  and  transmit  side,  even  when  considering  the  same  transmitter. 
Much  work  has  been  done  in  developing  more  accurate  models  for 
real  systems  (e.g.  [10]).  We  address  the  problem  of  general  corre¬ 
lation  structures  on  every  link  in  [8],  where  no  specific  correlation 
model  is  assumed.  A  jointly  optimal  solution  is  difficult  to  achieve, 
but  the  optimal  power  control  algorithm  presented  here  can  be  ex¬ 
tended  to  the  general  case,  and  a  very  good  suboptimal  solution  for 
the  transmit  and  receive  beamformers  is  also  given.  In  the  general 
case,  the  beamformers  are  able  to  intelligently  cancel  interference, 
and  we  see  that  savings  can  be  achieved  in  power  when  compared  to 
CSI  algorithms. 
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